Ionospheric TEC prediction using hybrid method based on ensemble empirical mode decomposition (EEMD) and long short-term memory (LSTM) deep learning model over India

Total electron data (TEC) from GPS nowadays can be used as a tool for understanding the space weather phenomena. The development of prediction model for TEC is quiet crucial and challenging due to the dynamic behavior of the ionosphere, since it depends on different factors such as seasonal, diurnal and spatial variations, solar geomagnetic conditions etc. In this paper, an attempt is made for predicting the GPS derived TEC values for different GNSS stations over India using a hybrid method based on Ensemble empirical mode decomposition (EEMD) and Long Short-Term Memory (LSTM) deep learning method. The daily TEC time series data from the IISc Bangalore (Latitude 13.021, Longitude 77.570), Lucknow (Latitude 26.912, Longitude 80.956) and Hyderabad (Latitude 17.417, Longitude 78.551) stations over India during the period 2008 to 2015 of solar cycle 23 and 24 is used for analysis. The assessment of model performance for testing predicted output compared with LSTM and EMD-LSTM models, and their comparison results show that the hybrid EEMD-LSTM model presents better than the other models.     

Reachable set estimation for spacecraft relative motion based on bang-bang principle

For spacecraft formation flight, the information of relative motion reachable set is very important, which can be used to predict the operating boundary of adjacent spacecraft and thus to ensure the safety of spacecraft operation. In this paper, we aim at developing a numerical method to approximate the reachable set for spacecraft relative motion. In particular, we focus on the quality of the approximation and the computational cost. Based on the bang-bang control principle, a polyhedral approximation algorithm is proposed to compute the reachable set of a relative motion spacecraft system. An inner approximation and an outer approximation of the reachable set for the system can be obtained. We prove that the approximation quality measured in Hausdorff distance can be guaranteed. The method is easy to implement and has low computational cost. Finally, the effectiveness of the algorithm is demonstrated by experimental simulation.     

飞机载体的杆臂效应对GPS测速精度的影响

杆臂效应是由于测量体的安装位置与运动载体质心不重合而引起的。从理论上分析了杆臂效应产生的原理,推导得出杆臂效应引起的速度误差公式。通过仿真试验和飞行试验数据分析了飞机载体不同运动状态下,杆臂效应对GPS速度测量精度的影响。通过对杆臂误差补偿后,GPS测速精度得到很大提高。在惯导系统飞行试验中引入该方法,可以科学合理地给出系统试飞评定结果。     

Precise orbit determination and precision comparison for FY3C and FY3D with receiver antenna GPS and BDS PCV using spaceborne BDS and GPS observation data

The FY3C and FY3D satellites were equipped with global navigation satellite occultation detector (GNOS) receivers that received both GPS and BDS-2 signals. For further improving precise orbit determination (POD) precisions, we estimated receiver GPS and BDS signal phase center variations (PCV) models with 2° and 5° resolutions and set the different weights for GPS and BDS-2 observations in the combined POD. The BDS-based POD precision using BDS-2 satellite antenna phase center offset (PCO) values from the China Satellite Navigation Office (CSNO) are not as accurate as those obtained from the International GNSS Service (IGS) Multi-GNSS experiments project (MGEX). The estimated receiver GPS and BDS PCV models with 2° and 5° resolutions were estimated from the GPS phase residuals of GPS-based POD and BDS phase residuals of combined POD, respectively. In most cases, the POD precisions using the estimated PCVs with 2° resolution are superior to those with 5° resolution. The precisions of the BDS-based POD and combined POD were both improved by introducing the receiver BDS PCV models. The weighting for GPS and BDS-2 observations can further improve the precision of the combined POD. The tested results of selected weights are better than those with equal weight in the combined POD. The experiment results show that orbital precisions of FY3C are worse than those of FY3D.     

Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

捷联式重力仪在海洋测量中的应用与数据处理

针对海洋重力测量对重力测量系统的要求,给出了一种捷联式海洋重力仪SAG-2M,论述了系统组成、工作原理、特点和数据处理流程。利用近期获得的某海域重力测量数据,评估了捷联式重力仪精度。测量结果表明,重复测线的重力异常曲线吻合度较高,交叉点不符值精度约为1.2mGal,满足海洋重力测量的指标要求。     

微动垫夹持刚度对微动疲劳寿命的影响分析

以单卡头式微动疲劳试验装置为对象，设计了不同微动垫夹持型式的微动疲劳试验装置，建立了两种试验装置的物理模型，结合理论分析和有限元数值仿真，从微动垫夹持刚度的角度分析了影响微动疲劳寿命的主要因素，并进行了试验对比验证。结果表明:夹持刚度的减小对SWT（Smith Watson Topper）微动疲劳损伤参量的最大值及位置影响较小，且使其在接触区域的分布趋于对称，可显著地降低相对滑移幅值，进而减小微动疲劳过程中的磨损程度；与夹持刚度较大的试验装置相比，夹持刚度较小的试验装置测试得到的微动疲劳试验寿命较长。      

异构多星编队的脉冲构型保持控制方法

近地轨道的双星编队通常设计具有自稳定性的编队构型参数初值，通过保持编队构型参数形成长期稳定的相对周期运动。针对编队中卫星数量增多产生的相对运动耦合问题，提出了基于Hill坐标和三角函数公式的多星相对运动分析方法。基于SAR载荷测量基线定义，结合多星编队构型参数的相对运动特性，提出了编队构型参数的设计方法，能够实现多星编队的最大有效基线组合。通过分析J2项摄动和大气阻力摄动的长期影响，研究了异构多星编队的相对运动衍化规律，提出了主从形式的脉冲偏置控制，能够有效保持针对异构多星编队设计的编队构型。通过面质比异构的四星编队控制仿真，验证了脉冲偏置控制形式下异构多星编队构型保持控制方法的有效性。     

基于单目视觉的非合作目标相对位姿测量

为实现对空间姿态翻滚航天器的在轨服务与维护以及对空间碎片的清理，需对其进行精确的相对位姿测量。针对相对位姿测量问题，提出了基于单目视觉与卡尔曼滤波的相对位姿测量方法。通过对特征点匹配算法进行调查，采用了具有尺度不变性与旋转不变性的尺度不变特征变换算法(SIFT)和加速稳健特征算法(SURF)的特征点提取方法，并对二者进行了对比，得到了二者分别适用的工况条件。通过对Kalman滤波算法进行研究，引入了相机偏置矩阵，设计了Kalman滤波器，解决了单目相机的距离模糊问题，估计得到了非合作目标的相对位姿、主惯量比以及特征点位置信息。经过仿真，姿态角度估计误差在稳定后低于0.3°，相对位置估计误差在稳定后低于0.5m，相较于真值，误差小于1.67%，主惯量比估计误差在稳定后低于0.01，特征点位置误差在稳定后低于0.005m。在引入相机偏置条件后，滤波状态变量均收敛，并得到具有足够精度的估计，成功解决了单目相机深度信息缺失问题。